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Control of the magnetic phase coexistence in NdFe 3 (B0 3 ) 4 by electric fields 
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We present a resonant x-ray diffraction study of the magnetic order in NdFe 3 (B03)4 and its 
coupling to applied electric fields. Our high-resolution measurements reveal a coexistence of two 
different magnetic phases, which can be triggered effectively by external electric fields. More in 
detail, the volume fraction of the collinear magnetic phase is found to strongly increase at the 
expense of helically ordered regions when an electric field is applied. These results confirm that the 
collinear magnetic phase is responsible for the ferroelectric polarization of NdFe3(B0 3 )4 and, more 
importantly, demonstrate that magnetic phase coexistence provides an alternative route towards 
materials with a strong magnetoelectric response. 

PACS numbers: 75.85.+t, 77.80.-e, 64.70.Rh, 61.05. 

Frustrated magnetic materials exhibit a large variety 
of intriguing physical phenomena, including the concomi¬ 
tant appearance of coupled ferroelectric and magnetic or¬ 
ders mm. In fact, interlinked magnetic and ferroelectric 
orders are very rare in nature. Their discovery in variety 
of frustrated magnets was therefore a surprise and gen¬ 
erated a lot of excitement ms], not only because these 
phenomena are very interesting from the viewpoint of 
basic research. There is also a significant technological 
potential Eng. Especially since a strong magnetoelectric 
coupling enables to store information in an extremely en¬ 
ergy efficient way mi. Indeed, earlier experiments could 
demonstrate modifications of a magnetic order related 
to the reversal of the electric polarization |MTT]. These 
studies, however, were concerned with electric field in¬ 
duced changes within a single magnetic phase or the ma¬ 
nipulation of phase stability close to a magnetic transi¬ 
tion [12]. 

In the present study we consider a different situation, 
where frustrated magnetic interactions cause two distinct 
ordered states to be metastable well below the observed 
critical temperatures, even in the absence of any first or¬ 
der transition. As a specific example, we investigated 
the electric field dependence of the magnetic order in 
NdFe 3 (B0 3 )4, which is a frustrated magnetic material 
with a large magnetoelectric response [13] • We show that 
a small perturbation created by an applied electric field 
allows controlling the stability of the coexisting collinear 
and helical magnetic orders. This implies that phase co¬ 
existence driven by magnetic frustration underlies the 
magnetoelectric response of NdFe 3 (B0 3 )4, providing an 
alternative route towards large magnetoelectric effects. 


The rhombohedral lattice structure of this material is 
illustrated in Fig.|T| where, for the sake of simplicity, only 
the two magnetic sublattices of Fe and Nd are shown. 
These two sublattices are magnetically coupled and un¬ 
dergo two magnetic phase transitions as a function of 
temperature USE!: upon cooling, commensurate mag¬ 
netic (CM) order sets in first at In ~ 30 K. In this phase 
the spins are ordered in a collinear fashion, forming fer¬ 
romagnetic (FM) a£>-planes, which are coupled antiferro- 
magnetically (AFM) along the c-direction. The magnetic 



FIG. 1. (Color online) Crystal structure of NdFe 3 (B0 3 )4. 
The material can be described using the rhombohedral space 
group R32 with hexagonal unit cell parameters a « 9.5 A 
and c ss 7.5 A at room temperature [14]. Only the magnetic 
sublattices of Fe (red) and Nd (purple) are shown. Fe is coor¬ 
dinated by oxygen octahedra, which form edge sharing helical 
chains along the c-direction. These chains and NdOe trian¬ 
gular prisms are connected in the ab-plane by BO3 triangles 
(not shown). 
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FIG. 2. (Color online) L-scans around the (0 0 4.5) position at 5.8 K, fitted with three Lorentzian squared peaks and a linear 
background (continuous lines). The the sum of the fitted curves agrees well with the experimental data (dots). The process 
induced by the L-field of E a = 8.8 • 10 4 V/m is reversible. 


modulation vector is commensurate, causing magnetic 
superlattice peaks at qcM= (0, 0, 3?r+3/2) with integer n. 
Below Tic ~ 15 K, a continuous transition into an incom¬ 
mensurate magnetic (ICM) phase occurs, where the FM 
moments of the a 6 -planes form helices propagating along 
the c-direction. The helical order is signaled by magnetic 
superlattice peaks at qicM= (0,0,3n + 3/2 ± S), where 
S increases continuously with further cooling due to a 
dramatically decreasing period of the spin helix ranging 
from 523 x c at 14 K down to 146 x c at 1.6 K. 

Interestingly, applying an external magnetic field B a ss 
1.3 T at 4.8 K along the a-direction induces an electric po¬ 
larization V a of up to 400 /rC/m 2 along the same axis pH] , 
This strong magnetoelectric coupling sparked signifi¬ 
cant interest in NdFe 3 ( 603)4 and immediately raised 
the question about its microscopic origin. Recent x-ray 
diffraction experiments already provided microscopic in¬ 
sight and implied that the finite V a is related to the 
magnetic field-induced CM order mm- In principle, 
the strong magnetoelectric coupling discovered in previ¬ 
ous studies should also enable to alter the magnetic order 
of NdFe 3 (BC> 3)4 by applying external electric fields. Yet 
no experimental studies of this effect have been avail¬ 
able up to now. We therefore performed high-resolution 
magnetic Resonant X-ray Diffraction (RXD) in applied 
electric fields in order to close this gap. 

The magnetic RXD has been performed at the res¬ 
onant scattering beamline of the European Synchrotron 
Radiation Facility (ESRF) in Grenoble, France (formerly 
installed at ID20) [T9]. A liquid helium Orange cryostat 
equipped with an E-field stick in-house built by one of 
us (C.M.)- was mounted on a six circle diffractometer 
as described in Refs. m and 1201 using an horizontal scat¬ 
tering geometry. The sample with a thickness of 1.3 mm 


along the a-direction was mounted between two vertical 
electrodes, leaving an additional gap of 0.5 mm between 
the top electrode and the NdFe 3 (B 0 3 )4 crystal. By ap¬ 
plying a voltage across the two electrodes, electric fields 
up to E = 8.8 • 10 4 V/m parallel to the crystallographic 
a axis were generated. We set the photon energy to the 
NdL 2 edge at 6.726 keV in order to probe directly the 
magnetic order of the Nd-sublattice, which, in turn, also 
reflects the magnetic order of the Fe-sublattice [T7]. For 
the L-scans shown in the following, the incoming light 
was horizontally ( 7 r) polarized and no polarization ana¬ 
lyzer was used. 

In Fig. [2] L-scans through the magnetic superlattice 
reflections measured at 5.8 K are presented. Prior to 
these measurements the state of the sample was prepared 
by cooling down to 5.8 K without electric field. The L- 
scan obtained for this zero field cooled state is shown in 
Fig.[2](a). In agreement with previous results both the 
(0, 0, 4.5) of the CM phase as well as the ICM-satellites 
at (0,0,4.5±5) are observed, revealing the coexistence of 
the CM minority and the ICM majority phase under the 
present conditions. It is likely that the observed residual 
intensity at (0, 0, 4.5) is due to pinned CM-domains that 
survive deep inside the ICM phase. This observation al¬ 
ready indicates that both phases have very similar free 
energies and that small perturbations suffice to drive the 
magnetic system from one phase to the other. 

The effect of an electric field on this phase coexistence 
can be observed in Fig.[2](b). Applying E a =8.8-10 4 V/m 
results in a transfer of intensity from the ICM satellites to 
the central CM-reflection. The application of the electric 
field therefore stabilizes the CM-phase and increases its 
volume fraction at the expense of the ICM phase. This 
conclusion is verified by the fact, that the sum of the scat- 
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tered intensities from the IC and the CM domains stays 
constant within 5%. However, the transformation from 
the ICM to the CM phase is not complete within the stud¬ 
ied E-field range, as incommensurate reflections retain a 
significant intensity. After the .E-field is switched off, the 
zero field state is recovered, as demonstrated in Fig. 1(c). 
Hence the effect of the applied E-field is reversible and 
our results do not show any significant hysteresis. 

The E-field also has a corresponding effect on the 
full width at half maximum (FWHM) of the superlat¬ 
tice reflections, which measures the average correlation 
length of the CM and ICM ordered regions. In detail, 
the experimental resolution in reciprocal space, as deter¬ 
mined using the neighboring (0,0,6) Bragg reflection, is 
0.001 x 27r/c, which is equivalent to a real space resolution 
of about 760 nrn. At 5.8 K and E a = OV/m the width 
of the ICM-satcllites is found to be resolution limited 
(FWHM=0.001 x 27r/c), implying a correlation length 
iiGM > 760 nm. Increasing the electric field to E a = 
8.8 • 10 4 V/m causes a clear reduction of £/cm to about 
476nm (FWHM=0.0016x27r/c). Conversely, the correla¬ 
tion length of the CM-ordered regions increases with E a 
from (cjj ~ 205 nm at E a = 0 V/m to (cm — 450 nm at 
E a = 8.8-10 4 V/m (FWHM decreases from 0.0037 x 27r/c 
to 0.0017 x 27r/c). The comparable values for (cm and 
£icm again show that phase transformation from ICM 
to CM is not complete at our highest electric field. 

Note also that the observed changes as a function of 
E cannot be caused by a heating of the sample. Heating 
would result in a strong change of S , which, however, is 
not observed. The changes in intensity and FWHM are 
therefore truly related to a change in the volume fractions 
occupied by the CM and ICM domains. The stabilization 
of the CM by an applied electric field further implies 
that the electric moment of this phase must be larger 
than that of the ICM phase. In fact, the comparison to 
macroscopic measurements indicates that only CM-phase 
is ferroelectric while the ICM phase does not posses a 
significant V [2T| . 

From the present experiments it is clear that an ap¬ 
plied E-field has a large effect on the magnetic order of 
NdFe 3 (BC> 3 ) 4 . The effects observed here indeed go be¬ 
yond earlier observations [Mj in that the population of 
two distinct, but coexisting phases is altered by an elec¬ 
tric field even far away from a magnetic transition. In¬ 
terestingly, the transition between the ICM and the CM 
phases as a function of E at constant temperature does 
not involve the continuous rotation of spins from the heli¬ 
cal towards the collinear alignment (5 remains constant). 
This is in stark contrast to the behavior which is observed 
as a function of temperature in zero applied fields. 

When applying a magnetic field at a fixed tempera¬ 
ture, we also find that 5 remains constant, i.e., in this re¬ 
spect the effect of electric and magnetic fields is the same. 
Nonetheless, this observation is surprising and sets the 
present results apart from previous work: as we demon- 



FIG. 3. (Color online) Comparison of experimental (sym¬ 
bols) and simulated (lines) Stokes scans for circular right (c r ) 
and circular left (ci) polarized incoming photons. The mea¬ 
sured superlattice peaks are (0 0 7.5) and (0 0 7.5+5) for CM 
and ICM, respectively. Temperature and fields for the shown 
measurements are: CM-T: 20 K, no field. ICM: 5.8 K, no field. 
CM-B: 2K, 2T. CM-E: 5.8K, 8.8- 10 4 V/m. 

strated earlier El. upon field-cooling to a fixed T<T)c, 
the value of 5 indeed does change with the magnetic field 
strength. This shows that the magnetic state at a given 
(B,T)-point in the phase diagram, depends on the path 
to this this point. 

In order to get more microscopic information about 
the magnetic state under the various studied conditions, 
we performed a polarization analysis for circular right 
(c r ) and left (ci) polarized incoming photons by means 
of so-called Stokes scans: for these scans, the polar¬ 
ization of the incoming and the detected photons are 
both controlled. The incoming beam polarization was 
set to c r or c\ using a diamond phase plate, whereas 
the polarization of the scattered beam was monitored 
using a Cu220-analyzer. The polarization state of the 
incoming radiation, as characterized by the polarization 
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parameters (Pi. Po. |Pq|) [22]. was (0.030,-0.006,1.000) 
and (0.042,-0.008,0.999) for circular left and circular 
right, respectively, showing the high degree of circular 
polarization. During a Stokes scan, the polarization of 
the scattered beam in the plane perpendicular to its di¬ 
rection is measured as a function of the angle 77, where 
77 = 0° corresponds to vertical and 77 = 90° to horizontal 
polarization. We integrated rocking scans of the analyzer 
crystal at each 77, to obtain the integrated intensity. Ex¬ 
amples and further details of the polarization analysis in 
magnetic RXD can be found in, e. g., Refs. l2^UMl 

In this way, we investigated four different states: (i) 
the CM phase in zero fields at 20 K (CM-T), (ii) zero 
field ICM at 5.8K (ICM), (iii) the CM-phase induced 
by applied magnetic field (CM-B) at 2 K and B a = 2 T, 
as well as (iv) the CM-phase induced by applied electric 
field (CM-E). 

The experimental data shown in Fig. [3] was then com¬ 
pared to model calculations, in order to extract infor¬ 
mation about the microscopic magnetic configuration of 
the sample. For modelling the magnetic RXD at the 
AFM superlattice positions we used the leading first or¬ 
der term of the spherical harmonic expansion of the scat¬ 
tering tensor, F W [23 ESI ■ Since the magnetic moments 
are parallel to the afeplane HM7] and there is a 3-fold 
axis along c, only one fundamental spectrum of F ^ en¬ 
ters, allowing us to apply the usual formalism given in 
Ref.27. Specifically, we use F^ oc e^-feTO*, , with ey*, and 
TO*, being the totally antisymmetric 3 rd rank tensor and 
the magnetic moment direction, respectively. Employing 
these magnetic scattering tensors and the experimental 
geometry, the structure factor tensors for the different 
spin-configurations were calculated and compared to ex¬ 
periment. 

For the CM-T phase we find that the polarization de¬ 
pendent magnetic RXD cannot be described by a mono¬ 
domain state. Instead, 3 domains are necessary to repro¬ 
duce the experimental data (top row in Fig. [ 3 ]). While 
these measurements do not allow to determine the ab¬ 
solute spin directions of the different magnetic domains, 
we find that they must be rotated by 120° about c with 
respect to one another. This is expected because of the 3- 
fold axis of the unperturbed lattice structure (neglecting 
the small ferroelectric distortion). 

The onset of the ICM order is signaled by the appear¬ 
ance of the IC satellites, which exhibit different resonant 
intensities for c r and c\ (second panel from top in Fig. [3]). 
This difference is due to the chirality of the magnetic he¬ 
lix. Specifically, our analysis yields the expected presence 
of left and right handed domains with volume fractions 
of 0.70 and 0.30, respectively. 

The above findings for the CM-T and ICM phase are in 
good agreement with neutron scattering results [T5i in. 
It was also shown earlier that a magnetic field of 2T 
applied to the CM phase suffices to create a mono-domain 


CM-B state [171 IB] ■ Our analysis of the Stokes scans 
yields the same result, using a different and independent 
method. The above comparison to previous studies hence 
verifies our approach and shows that the full polarization 
control in magnetic RXD provides a powerful means to 
identify and characterize magnetic multi-domain states. 

Turning to the CM-E, we find that 3 magnetic domains 
with a relative angle of 120° are necessary to model the 
experimental data, very similar to the case of the CM-T 
state described above. We therefore find that our max¬ 
imum electric field of 8.8 • 10 4 V/m does not create a 
mono-domain state with collinear magnetic order along 
one single direction. Instead still different types of do¬ 
mains are found to exist. This is in contrast to the mono¬ 
domain CM-B state described above. 

In summary, the key result of the present work is the 
observation of a delicate balance between two distinct 
magnetic phases, which reacts very sensitively to an ap¬ 
plied electric field. More in detail, the present magnetic 
RXD experiments reveal a complex microscopic state, 
where the coexisting CM- and ICM-regions also realize 
different domain-types, namely three 120° CM- and two 
left/right handed ICM-domains. It is important to point 
out that frustrated magnetic couplings result in the com¬ 
petition of the CM and ICM order D33- Hence, this frus¬ 
tration also underlies the magnetic phase coexistence re¬ 
ported here, which occurs even deep inside the ICM phase 
and in the absence of a first order transition. Triggering 
the magnetic phase coexistence in frustrated materials 
by electric fields therefore provides an alternative route 
towards new materials with a large magnetoelectric re¬ 
sponse. 
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